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Abstract 

We study a discrete model of the Laplacian in R 2 that preserves the 
geometric structure of the original continual object. This means that, 
speaking of a discrete model, we do not mean just the direct replacement 
of differential operators by difference ones but also a discrete analog of 
the Riemannian structure. We consider this structure on the appropriate 
combinatorial analog of differential forms. Self-adjointness and boundness 
for a discrete Laplacian are proved. We define the Green function for this 
operator and also derive an explicit formula of the one. 

1 Introduction 

We begin with a brief review of some well known definitions that are related 
to the contents of this paper. Denote by A r (IR 2 ) the set of all differentiable 
complex-valued r-forms on M 2 , where r = 0, 1, 2. Let * : A r (R 2 ) -> A 2 ~ r (]R 2 ) be 
the Hodge star operator. An inner product for r-forms with compact support 
is defined by 



O, ip) = J pA*ip, (LI) 

R 2 

where the bar over if; denotes complex conjugation. Let L 2 A r (IR 2 ) denote the 
completion of A r (R 2 ) with respect to the norm generated by the inner product 
(fTTTj) . Let the exterior derivative d : A r (R 2 ) -> A r+1 (M 2 ) be defined as usual. 
We define the operator 

d : L 2 A r (R 2 ) -> L 2 A r+1 (M 2 ) (1.2) 
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as the closure in the L 2 -norm the corresponding operation specified on smooth 
forms. The adjoint of d, denoted by 5, is given by 

(dip, w) = (p, 6lj), ipeL 2 A r (R 2 ), lu e L 2 A r+1 (R 2 ). 

Note that 6 : L 2 A r+1 (M. 2 ) -> L 2 A r (R 2 ). The following relations hold among *, 
d and 6. See for instance [T5] . 

^2 = (_!)r(2-r) 7d) 5 = (-l) r t' 1 d * . (1.3) 

The Laplacian is defined to be 

- A = dS + Sd : L 2 A r (R 2 ) -> L 2 A r (R 2 ). (1.4) 

In this paper we develop some combinatorial structures that are analogs of 
objects in differential geometry. We are interested in finding of a natural discrete 
analog of the Laplacian on cochains. Speaking of a discrete model, we mean not 
only the direct replacement of differential operators by difference ones but also 
discrete analogs of all essential ingredients of the Ricmannian structure over a 
properly introduced combinatorial object. 

Our approach bases on the formalism proposed by Dezin [5]. We adapt 
the combinatorial constructions from [HJ [12] and define discrete analogs of 
operators (ll.2| - (ll.4|) in a similar way. In [5], Dezin study discrete Laplace 
operators in finite-dimensional Hilbert spaces, i.e. on cochains given in domains 
with boundary. In this paper we extend these results on an infinite complex 
of complex-valued cochains. We prove self-adjointness and boundness for the 
discrete Laplace-Beltrami operator in infinite Hilbert spaces that are associated 
with L 2 A r (R 2 ). Spectral properties are discussed. We define the Green function 
for the discrete Laplacian and derive the one in an explicit form. 

There are other geometric approaches to discretisation of the Hodge theory 
of harmonic forms presented in [6], [7l [8] . In all these papers discrete models are 
given on the simplicial cochains of triangulated closed Riemannian manifolds. 
See also [21 El [9] and references given there. 

Classical references on second order difference equations are the books by 
Berezanski [2J ch. 7], Atkinson [Tj and the most recent monograph by Teschl 

2 Combinatorial structures 

Let us denote by £(2) the two-dimensional complex. This complex is defined by 

c(2) = £° © e 1 © £ 2 , 

where £ r is a real linear space of r-dimensional chains. We follow the notation of 
[51H1]- Let {xk, s }, i e l s > e fcsi' {^fc,s}, fc, s € Z, be the sets of basis elements 
of £ , £ , £ 2 respectively. It is convenient to introduce shift operators 

rk = k + 1, ak = k — 1 



2 



in the set of indices. The boundary operator d is defined by the rule 

9x k , s = 0, dfl ktS = e\ s + e 2 Tk s - e\ TS - e 2 ks , 

del s = x rk,s - Xk,s, del >s = x k ,Ts - Xh,s- (2- 1 ) 

The definition of d is linearly extended to arbitrary chains. We call the complex 
£(2) a combinatorial model of R 2 . 

On the other hand, we can consider £(2) as the tensor product £(2) = £© £ 
of the one-dimensional complex £ (combinatorial model of a real line). Then 
basis elements of £(2) can be written as follows 

x k <£> x s = x k<s , e k ® x s = el s, 
e k © e s = tt k:S , x k (g> e s = e 2 ks , 

where x k , e k are the basis elements of £. 

Let us introduce an object dual to £(2). Namely, the complex of complex- 
valued functions over £(2). The dual complex K(2) we can consider as the set 
of complex- valued cochains and it has the same structure as £(2), i.e. K(2) = 
K° © K 1 © K 2 . In other words, K{2) is a linear complex space with basis 
elements 

{x k > s , e k ' s , e\ s , n k > s }. 
The pairing (chain-cochain) operation is defined by the rule: 

< x k , s , x™ >=< Cl k>8 , n p ' q >=< e 1 ^, e v { q >=< e\ s , e p 2 q >= 5™, (2.2) 

where 5 k ' q s is Kronecker symbol. We call elements of the complex K(2) forms, 
emphasizing their closeness to the corresponding continual objects, differential 
forms. Then the 0-, 1-, 2-forms (p, u) = (u, v), rj can be written as 

V> = X! L Pk, s x k ' s , rj = ^ri kyS n k ' s , lu = ^(u k!S e'l' s + v k;S e 2 ' s ), (2.3) 

where tp k ,s, u kySl v k>s , f] k)S € C for any fc, s G Z. Operation (|2.2[) is extended to 
arbitrary forms (|2.3[) by linearity. The boundary operator (|2.ip in £(2) induces 
the dual operation d c in K(2): 

< da, a >=< a, d c a >, (2.4) 

where a S £(2), a £ if (2). We assume that the coboundary operator d c : K r — > 
X r+1 is a discrete analog of the exterior differentiation operator d (|1.3p . 

If y> <G X° and uj = (u,v) E K , then we have the following difference 
representations for d c : 

< el >s , d c (p >= (p TkiS - tp k<s = A k tp k , s , 

< fl k , s , d c u >= v TktS - v kiS - u ktTS + u kiS = A k v k , s - A s u k:S . (2.5) 
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Note that if 77 e K 2 , then d c r] = 0. 

Let us now introduce in K(2) a multiplication which is an analog of the 
exterior multiplication A for differential forms. We denote this operation by U 
and define it according to the rule: 



x k ' s 


U x k ' s = 


k.s k,s 

x , e 2 


1 1 k,rs nk.s 

U e{ — —il ■ , 




x k ' s 


U e k { s = 


e k,s y x rk,s = 


k.s 

ei , 




x k > s 


U e^ s = 


e k ' s U x Kts = 


k,s 
e 2 -> 




x k ' s 


U n k ' s = 


: n k ' s ux rKrs 


= e k ' s U e T 2 k ' s = fl k ' s , 


(2.6) 



supposing the product to be zero in all other cases. The U-multiplication is 
extended to discrete forms by linearity. In terms of the theory of homologies, 
this is the so-called Whitney multiplication. For arbitrary forms a, j3 S K(2) 
we have the following relation 

d c (aU/3) = d c aU{3+ (~l) r aUd c f3, (2.7) 

where r is the degree of a. The proof of this can be found in Dezin p. 147]. 
Relation (|2.7p is an analog of the corresponding continual relation for differential 
forms (see [15]). 

Define a discrete analog of the Hodge star operator. Let e k,s denote an 
arbitrary basis element of K(2). We introduce the operation * : K r — * K 2 ~ r by 
setting 

e k < s u*e M = n k ' s . (2.8) 

Using (|2.6[) we get 

*x k ' s =n k ' s , *e\ s = e rk '\ * e k 2 ' s = -e k { TS , *n k ' s = x Tk > TS . 

The operation * is extended to arbitrary forms by linearity. 
Let a € K r is an arbitrary r-form: 

a = Y,Oik, s e k ' s . (2.9) 

k,s 

Denote by Kq the set of all discrete r-form with compact support on £(2). Let 
fi be the following "domain" 

fl = J2 Q k, s , k,seZ, (2.10) 

where Qk, s is a two-dimensional basis element of £(2). Note that if the sum 
(|2~TU)1 is finite and let -N < k,s < N, N € N, then we will write Q = Q N . 
The relation 

(a, (3) =< fi, qU^>, (2.11) 
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where a,/3£ Aq, gives a correct definition of inner product in A(2). Using 
(|2.2p . (|2.6p and (|2.8[) , this definition can be rewritten as follows 

(a, j9) = 5^a M ^7. (2.12) 

For = O^v we will write 

jv 

(a, (3) N =< Qjv, a U */3 >= ^ a ktS (5 k ^ s . 

k,a=-N 

Let a e if r , /? £ A r+1 . The relation 

{d c a, P) N =<dn Nl aU^>+(«, <5 c /3)at (2.13) 

defines the operator 5 C : A r+1 —> A r , 

£ c /3 = (-l) r *" x d c * /3, 

which is the formally adjoint operator to d c (see [5] for more details). It is 
obvious that the operator 5 C can be regarded as a discrete analog of the cod- 
ifferential 8 (cf. (II. 3|) ). Equation (|2.13|) is an analog of the Green formula for 
the formally adjoint differential operators d and 5. It is easy to check that for 
a G A„, (3 S Kq +1 we obtain 

(cf a, /3) = (a, (5 c /3). (2.14) 

According to (|2.5p . we have i5 c (/? = and 

< x fc , s , 6 c uj >= -A k u ak ,s - A s v k ^ s , (2-15) 

< e£ iS) <S c r/ >= A s r)k,a S , < e\, s , ^ >= -A k r] aktS , (2.16) 

where y> € A , uj € A 1 and 77 e A 2 . 

Therefore a discrete analog of the Laplace-Beltrami operator (|1.4p can be 
dchncd as follows 

- A c = 5 c d c + d c S c : K r -» IT. (2.17) 
Obviously, if <p € A" , then we have 

-6?<p = 8 c dp<p. (2.18) 

Combining (|2.15p with (|2.5p we can rewrite (|2.18p as 

< Xfc )g , ~A C tp >= Aifk,s - ip Tk}S - ip k , TS - ^ CTfe)S - fk,as- ( 2 -19) 

The same difference form of (|2.17p can be drawn for the components rjk^ s of 
r) £ A 2 and for the two components Uk, s Vk, s of w € A 1 . 
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3 Discrete Laplacian 

Let us now introduce the linear space 

H r = {aeK r : ^K s | 2 <+oo, k,seZ}, (3.1) 

k,s 

where r = 0, 1, 2. Clearly, TC r is a Hilbert space with inner product (|2.1ip (or 
(|2.12p ) and with the following norm 

\\a\\ = ^a) = (Y.M 2 Y ■ ( 3 - 2 ) 

k.s 

Note that if a £ 7i r , then the set of complex- valued sequences {ctk, s ) is l 2 (l?). 
From now on we regard d c , S c and — A c as the following operators 

d c :H r ^H r+ \ 6 c :H r ^H r - 1 , -A c :H r ^H\ 

where r = 0, 1, 2. It is convenient to suppose that Or 1 = H? = 0. 

Theorem 1. The operators 

-A c :H r ^H r , r = 0,1,2, (3.3) 

are bounded and self-adjoint. Moreover, \\ — A c || = 8, where \\ — A c || denotes 
the operator norm of — A c . 

Proof. We begin by proving self-adjointness of — A c for the case r = 0. Let 
(p, ip £ Tt° and lu = (u, v) € 7i . Then the Green formula (|2.13|) can be rewritten 

as 

JV 



(d c cp, uj) n = {Vk,rNVk,N - (fk-NVk- T N) + 

k=-N 
N 

+ 22 ( 1 PtN,sUN^ ~ (p-N,sU-TN,s) + (<P, 5 c lo) n (3.4) 



s=-N 

The substitution of d c ip for u> in (|3.4p gives 

JV 



(d C lfi, d C 1p)N = 22 ( l fik,TN(lpk,TN - tpk,N) - ¥>k,-N(lpk,-N ~ ^fc,- T jv)) + 
k=-N 

N 

+ {VrN,s{lj)TN,s - 1pN,s) ~ <£-iV,s(V>-iV,s ~ Ip-rN,*)) + {<P, S C d C ^) N . 

— N 

(3.5) 



Letting N — > +oo we get 

<fty) = (<p, S c d c iP). (3.6) 
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It follows immediately that — A c : H° — > H° is self-adjoint. 

The same proof remains valid for the case r = 2. Now we have — A c = d c S c 
and the analog of relation (|3.5p results from the inner product (6 c rj, S c ()n, 
where jy,( € Tt 2 . A trivial verification shows that properties of (S c r), 6 c Qn are 
completely similar to those of (d c ip, d c V>)iv- Hence 

(5% SX) = (r), d c 5 c (). (3.7) 

Finally, let r = 1. In this case we have — A c = d c 5 c + S c d c and we must 
study the sum (d c uj, d c d) + (5 c uj, <5 c tf), where uj = (u,v) £ H 1 , & = (f,g)eH l . 
Taking in (|2 . 1 3[1 a = ui and fj = d c d we obtain the analog of relation (|3.5[) for 
1-forms 

N 

Tk,s) 



(d c u>, d c i!)) N = ^ (AfeW)t, s - A s Ufc iS )(Afcgfc, s - A a /; 

k,s=-N 

N 

= ^ [uk,-N{Ak9k,-TN — A-Nfk,-rN) — Ufc,TN(Ak9k,N — Afffk,N)] + 
k=~N 
N 

+ [ t, TJV, s (AjvS'jv, s - A s /at iS ) - v-N iS (A N g_ rNt3 - A s /_ tA t iS )] + 

s=-iV 

Letting TV — > +oo we obtain equation (|3.6|) for the 1-forms uj, $ £ H 1 . In the 
same manner we can see that equation (|3.7[> holds for w, -d £ H 1 . Adding we 
obtain 

{d c uj, d c d) + (5 c uj, (5 c tf) = (uj, -A c tf). (3.8) 

Thus it follows that 

(-A c uj, i?) = (uj, -A c $). 

For the rest of the proof let a £ TL r be an arbitrary r-form. Substituting 
(Prgj) into we get 



0I = | ^(4a/c,s - Brts - Otk.rs - Ot a k.s ~ Oik,as)a>k r 
k,s 

< 4^ \ak,s\ 2 + X! \ a Tk,sUk^\ + 2_j \ a k,T S Ok^\-\ 



< 



^ |a<rfc, s afc :S | + |afc,<r S afc :S | < 8||a|' 2 
k , s 



From this we conclude that || — A c || < 8. Since — A c is self-adjoint, it follows 
easily that || - A c || = 8 (see for instance 10, Ch. 3]). 

□ 

Corollary 2. The operators (3.3\) are positive, i.e. for any non-trivial r-form 
a £ TL r we have 

(-A c a, a) > 0. 
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Proof. This follows from (gU), (02]) and (03). □ 
Corollary 3. For any r, r = 0, 1, 2, we have 

a(-A c ) = [0, 8], 
where o~{— A c ) denotes the spectrum of — A c . 

Proof. Straightforward. □ 

4 Discrete analog of the Green function 

Let g(— A c ) = C \ a(— A c ) denotes the resolvent set of — A c . In this section 
we try to describe the resolvent operator (— A c — A) -1 , A € g{— A c ), of the 
operator — A c : 7i° — > 7i . Let us introduce a discrete form G(x,x) on £° x €° 
as follows 

G(x, x) = Gk, s (x)x M , where G M (x) = ^ G fe , s , m , n x ro ' n 

and Gk, s .rn,n S C for any fe, s, m, n G Z. Hence we have 

G(x, 5) = ^ £ G M , m ,„x fe < s x m '". (4.1) 

k,s m,n 

This is a so-called discrete double form (for details see [11]). Note that in 
the continual case (for differential forms) this construction is due to De Rham 
|13J . It is obvious that the basis elements x k,s and x m,n in (|4.ip commute and 
G(x,x) = G(x,x). By analogy with (|2.2[) . the double 0-form G(x,x) can be 
written pointwise as 

^ •Ek,S'E J m i m G(x, x) ^> Gfc ;S , m;n . (4.2) 

Let <p E H°. Then we define a 0-form 8 m ' n , m,n € Z, by setting 

(p, ^'") = X>mC" = (4.3) 

where (5^™ is the Kronecker delta. By analogy with the continual case, the 0- 
form 5 m ' n defined by (|4.3p will be called a discrete analog of Dirac's (^-function 
at the point x m<n . We can write 5 m,n as 

k,s 

We need also the following double form 

d(x, x) = 2^ 2^ 5 fe,« 25 21 = 2^ 5 x = l^ x x ■ ( AA > 

k,s m,n m,n m,n 



8 



Definition 4. The double form is called the Green function for the oper- 

ator -A c : H° -> H° if 

G k , s , m , n (X) = (<5^ s , (-A c - X)~ l 5 m,n ) (4.5) 

for any k, s,m,n £ 1. 

Of course, 

(-A c -A)- 1 ,5 m '"-G m ,„(x,A). 

It follows easily that 

(-A c - X) x G(x, x, A) = 5(x, x), (4.6) 

where (— A c — A)^ is the operator — A c — A that acts with respect to x. Recall 
that in our abbreviation x corresponds to k, s. 
We have 

(-A c - A)- V = J2 Gk mi (^ m ,nx k ' s x m ' n , veH°, A e g(-A c ). 

k,s m,n 

Indeed, applying (— A c — X) x to the right-hand side gives 
( — A c — A)a; ^ ' l Pm,nG m , n (x, X)x m ' n = S ' <^ m ,n( — A c — X) x G m , n [x, X)x m ' n = 

m,n m,n 
m.n m,n 

We now try to write the Green function for — A c in a somewhat more explicit 
way. For this we construct a solution of the equation 

- A c ip = Xip, A e C. (4.7) 

The following construction is adapted from [14], where the Green function is 
studied for Jacobi operators. By (|2.19p , equation (14. 7p can be written pointwise 
(at the point Xk, s ) as 

4<y9fc, s - (f T k,s - <Pk,TS - <fak,s - <£k,as = X<Pk,s- (4-8) 

Applying the transformation A = —4^ + 4 we reduce (|4.8p to the equation 

^{<Prk,a + fk,rs + <Pak,s + <Pk,os) = Wk,s- (4-9) 

An easy computation shows that, substituting the ansatz ipk,s — p k+s into (|4.9p . 
we obtain 

^ s (/x) = ( M ±fl(M)) fe+s , (440) 

where R{fx) = — \] — 1 and denotes the standard brunch of the square 
root. It follows that we can write the solutions >p (X) of (14. 7p in the form 

^(A) = 5>±>)z fc ' s , (4.11) 
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where tpf s (/j,) are given by (|4.10[) and ^ = 1 — j. Obviously, A e [0, 8] if and 
only if /i e [-1, 1]. By Corollary 3, A e g(-A c ) leads topeC\ [-1, 1]. 
It is convenient to write ^ s (m) as 

V*» = ^(M)-^(M), (4-12) 

where (pkit 1 ) = U* ^ R(l J ')) k an d k,s <E Z. Suppose /Lt G C \ [—1,1]. Then one 
has to examine that the sequences (tp k (fJ*))kez, are square summable near ±oo 
respectively, i.e these are £ 2 (Z) near ±oo. 

Let H± denotes the set of 0-forms whose restriction to K± belongs to TiP . 
Here K\ (K°_) denotes the set of 0-cochains (f23|) with k + s > (ft + s < 0). 

Proposition 5. Let A S g{~A c ). Then <^ ± (A) S W±. 

Proof. This follows immediately from (|4. 12[) . □ 

Lemma 6. For any k,s £ 1 the component f^(n) ■ tp~ (n) is a solution of 

Proof. It is easy to check that 

<pf - 2(upf + 1 = 0. 

It follows that <Pk{fj) ■ <pj(n) satisfies equation (|4.9p . Indeed, putting this in 
(|4.9p . we have 

□ 



Theorem 7. Let X 6 g{— A c ). Then the components !j4-5\ ) of the Green function 
are given by 



Gk,s,m,n (A) 



4i?(/i) 



(/i + i?(M)) |r ^ m|+|TS ^ l| / rfc = ™, s>n 

or k > m, s = n, 

( M + fl(M)) kfe " m|+k ^ n| fork = m, s<n (4.13) 

or k < m, s = n, 

(/i + R(n)) ' fc m ' + s 71 for the all others, 



where (i = 1 — ^ . 



Proof. We must prove that 

< z M z m ,„, (-A c - A^Gfo 5) >= £™f , (4.14) 

where G(x,x) is given by (|4.1[) . Using (|2. 19[) and (|4.ip , we can rewrite the 
left-hand side of (|4. 14|) as 

(4 A)G/ Cs _ mn G T fc srnn Gak ,s,m,n G/;;_ rS)m _ n G^(j S;mj)l 

— 4:{J<Gk,s,m,n (Grk^s^n.n ~\~ Gak,s,m,n ~t~ Gk,TS,m,n ~\~ Gk,as.7n,n ) ■ (4-15) 
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The proof falls naturally into three parts. 

Fix x m>n . First, let k = m and s = n. Substituting (|4. 13|) into (|4.15|) and 
using (|4.12|) , we obtain 

< Xk, s X m ,n, (~A C - X) x G(x,x) > = j^(-fj,<p£([ji)tp£ 0) + (pf(fl)lf^((l)) = 



= 1. 



Note that tpf(p) = (ji ± R(ji)) Q = 1. 

Now we show that (14. 13)) satisfies equation (|4.8[) for the cases k — m, n 
and fc ^ m, s — n. Check the case k — m, s > n. In this case the left-hand 
side of (|4.14[) is equal to 

- <ptto<Pi-nM ~ ¥>t(»)vt-n(»)] = - W + 1) = 0. 

The proof for the case fc = m, s < n (or fc 7^ to, s = n) is similar. 

Finely, let k ^ m, s ^ n. We give the proof only for the case fc > to, s < n. 
For the other cases the proof runs similarly. In this case we have 

G k , s ,m,nW = qirfLWrf-.W = _4R( M ) VkitiVm O^ft (tivT (m) • 

Here we use that = y>r for any fc G Z. By Lemma 6, (p^ (/i)ipj (fj,) is a 
solution of (|4.9|) and so is Crfc, s ,m,n(A) at the point Xfc iS . 

Thus, since G(x,x,\) with components given by (|4.13|) and, by Proposi- 
tion 5, G(x,x, A) G 7i° with respect to x, it must be the Green function of 
-A c : H° -> H°. □ 

It should be noted that the consideration of — A c : Ti 2 — > Ti 2 does not differ 
from that carried out for the 0-forms ip G H°. In this case we consider (|4.1[) as a 
double form on £ 2 x £ 2 with basis elements r2 fc,s fl p ' 9 . Then we define the Green 
function as above and its components are given by ()4.13j) . The situation with 
— A c : H 1 — > H 1 is more difficult. In this case, having written equation ()4.7|) 
pointwise (at the elements e\ s and e? s ), we obtain a pair of equations (|4.8[) 
which corresponds to the components Uk, s , Vk,s of u> G W 1 . Roughly speaking, 
here we must describe the Green function for each component of the 1-form 
lo = {u, v). The reader can verify that the components Gfc, s ,m,n(A) of the Green 
function have again the form (|4~13)) . 
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